Composite bosons in bilayer v — 1 system: An application of the Murthy-Shankar 

formalism 



Ivan Stanic and Milica V. Milovanovic 
Institute of Physics, P. O. Box 68, 11080 Beograd, Serbia and Montenegro 
(Dated: February 2, 2008) 

We calculate the dispersion of the out-of-phase mode characteristic for the bilayer v = 1 quantum 
Hall system applying the version of Chern-Simons theory of Murthy and Shankar that cures the 
unwanted bare electron mass dependence in the low-energy description of quantum Hall systems. 
The obtained value for the mode when d, distance between the layers, is zero is in a good agreement 
with the existing pseudospin picture of the system. For d nonzero but small we find that the 
mode is linearly dispersing and its velocity to a good approximation depends linearly on d. This 
is in agreement with the Hartree-Fock calculations of the pseudospin picture that predicts a linear 
dependance on d, and contrary to the naive Hartree predictions with dependence on the square-root 
of d. We set up a formalism that enables one to consider fluctuations around the found stationary 
point values. In addition we address the case of imbalanced layers in the Murthy-Shankar formalism. 



A major problem surfaced in the early Chern-Simons (CS) composite boson description yj of v = 1 quantum 
bilayer. Namely, a bare electron mass appeared in the Bogoliubov out-of-phase gapless mode dispersion, which is 
unwanted due to the expectation that the leading description of any intra-Landau-level collective mode depends solely 
on interactions. The same requirement applies for the description of any quasiparticles that may exist in the lowest 
Landau level, namely their mass should stem from interactions Q. Murthy and Shankar 0,Q put forward an extended 
CS theory that was able to prescribe such a mass for composite particles. The theory underwent major additions and 
improvements (including an extension to higher momenta) Q , but in this paper we will use the early version £| best 
suited for our needs, i.e., the calculation of the dispersion relation of the out-of-phase mode in a small-momentum, 
low-energy window. 

On the other hand, what is believed to be the first calculation of the dispersion of the gapless mode in the scope 
of the pseudospin picture (in which an electron can be in states that are superpositions of localized layer states) was 
done by Fertig 6J. He obtained an interaction-dependent velocity of a linearly dispersing mode when d is nonzero 
and a spin wave quadratic dispersion with interaction-dependent spin stiffness when d is zero. The problem at d = 
was again addressed in Ref. |7| where a quantum ferromagnet picture for this case was established. Then followed the 
pseudospin picture of reference Ref. a standard reference for the bilayer v = 1 problem, in which small- momentum 
dispersion of the gapless mode was in agreement with Fertig's. 

Our Hamiltonian is 

U = l^ 2^ + Ve + V a , (1) 

1,(7 

where A is the vector potential of the constant external magnetic field —B a e z , so that the average total density is 
n = v/2nl B with v = l/(2s + 1), s = 0, 1, . . . (for generality) and Ib = l/\/eBo, the magnetic length. We take Ib = 1 
and h = l. m is the "bare" electron mass, which is precisely the effective mass of electron in GaAs. Va and Vb denote 
intralayer and interlayer interactions, respectively, a =T,| is the layer index. 

What follows is a brief account of a simple generalization of the Murthy and Shankar CS theory extended to the 
case of the bilayer. Missing details and explanations that are relevant also to the single-layer case can be found in 
Rcfs. 0,0, 01- As usual in the CS theory we make the unitary transformation 

C/ = exp|z(2 S + l)^$ y |, (2) 

where $y is the phase of the difference, Zi — Zj, of any two coordinates zi and Zj. Therefore, we emphasize the 
transformation is the same, irrespective of the layer indices. The corresponding Hamiltonian is 

Has = U-^HU = V + + Scs]2 +V E + V A (3) 
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where the new "gauge" field satisfies the following connection with the total density, p(r): 

V x a CS = 2tt(2s + l)p(r)e z . (4) 
Then we consider averaged and fluctuating values of acs and p 0, rewriting Eq. Q as 

V x (a C s) + Vx : a C s ■= 2n(2s + l)n + 2tt(2s + 1) : p : (5) 
so that equivalently, due to the cancelation of the external and averaged CS field, we can rewrite Eq. as 

Has = Y d Wi "' + 'f CS ' A2 +V E + V A . (6) 



I, a 



Now Shankar and Murthy 0,0, analogously to what Bohm and Pines 01 did in the case of a three-dimensional (3D) 
Coulomb gas, introduce magnctoplasmon degrees of freedom as separate and elementary but necessarily satisfying 
certain constraints with particle degrees of freedom in order to avoid overcounting. They do this by introducing a 
pair of conjugate fields, a(q) and P(q), for each q in a disk in the momentum space, 

[a(q),P(q')]=i(2Tr)H(q + q'). (7) 

[reminding us (for fixed q) of the harmonic oscillator commutation relation in (x,p) representation] and further, 
defining a longitudinal and a transverse field, P(q) and a(q), respectively, as 

P(q) = iqP(q) and a(q) = -ie z x qa(q). (8) 
They rewrite the Hamiltonian (density) as 

H= ^J2 *csJ-^+ : a cs ■ +S) 2 ^cs,a + V A + V E , (9) 

a 

in the second-quantized language, with the requirement (constraint) that 

a(q)|physical state) = 0, (10) 

for each q such that \q\ < Q where Q is the radius of the disk. Now it is convenient to eliminate : acs'- in favor of a 
and P, and Shankar and Murthy do that 3] applying the following unitary transformation, 

U = exp [i £ P{-q) M2S + l) p{q\ (11) 
1 ife q J 

where p(q) = p-\(q) + Pi(q) i.e. the total charge density. Now for the Hamiltonian density we have 

H = Y. ^*fcMHV + a + 27r(2s + l)P + M)H kb , a + V A + V E (12) 

where "ffc^o-'s denote the transformed fields, which describe modified, transformed quasiparticles, and 5a is the remnant 
of the CS field left uncanceled for \q\ > Q. The constraint gets the following form: 



2tt(2 S + 1) . 
a{q) p{q) 

Neglecting Q 5a from the start and introducing Q 



| physical state) = for |<f| < Q. (13) 



= a(q)+i2n(2s+l)P( q l 
W ' y/4ir(2s + 1) V ' 

and 

c (q) = Q- ^2Pj,<?+ exp{-iq- f jt(T }, (15) 
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where V± = V x ± iV y for an arbitrary vector V = V x e x + V y e y , we can rewrite our Hamiltonian as 



n = E + E "^(Mfl + J < 2s + V X! ftfl^® + c + (g)A(5)] + V A + V E . (16) 
^ — ' Ira *■ — ' V m L — ' 

*,<> \q\<q m<Q 

As expected oj c = eBo/m, i.e., equal to the cyclotron frequency and in deriving the magnetoplasmon term we neglected 
also total density fluctuations |4|- 

Our quasiparticles are bosons and again for the sake of completeness and easy reference we give brief account of the 
so-called final representation in the Murthy-Shankar approach applied to the case of two species of composite bosons. 
First Murthy and Shankar always approximate as 

exp{i(q - k) ■ r* } w n (27r) 2 <S 2 (<f- k) (17) 

in the long-wavelength approximation so that also in this bosonic representation we consistently have 

[c(S),ct(?')]«0. (18) 

To decouple the oscillators and particles they apply the following canonical transformation, 

U(X ) = exp{*SoAo} = cxp{A o ]T [J%A% - AHq)c(q)}}, (19) 

l<zl«2 



where 



(20) 



^4^(25 + l)n 

and the parameter Ao should be determined. As we have new variables f2 defined through 

n old = exp{-iS , A }^exp{iS'oAo}, (21) 



Murthy and Shankar also define 



so that in this case we have 



0(A) = exp{-iS , oA}^exp{iS , A}, (22) 



A(q*,X)=A(q)-6Xc(q), 

c(q,X)=c(q), \q\<Q. (23) 

It is easy to check that Ao = 1 does the job of decoupling and now we concentrate how variables p° ld (q), a =\, j arc 
connected with new ones. We use their definitions p^ ld (q) — J^expjig- r*i irT }, a =|, J, apply Eqs. (|21|l and l|22|l . find 
dp a (q, X)/d\, and integrate over A to get 



P°J d (q) = P*(q) + 



[A{q) + A\-q)]--[c{q) + c\-q)] 



(24) 



y/4n(2s + 1) n 

Immediately we can conclude that the spin density p° ld (q) = p° ld {q) — P° ld (q) is invariant, 

pf\q) = pM, (25) 

under the transformation of the final representation [Eq. (|19fl ]. when we assume that we have the same fixed number 
of particles, = = n/2, in each layer. This is our main claim in the Murthy-Shankar formalism for the bilayer 
system. 

The analysis for the charge density and the form of the constraint in new variables proceeds as in Refs. 0,0] and 
finally for the form we get 

p(q) = -|E^ X exp{~iq- fj,„}. (26) 
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In a few lines but also using an assumption that we deal with an infinite system (without boundary) we can prove 
that, in the second-quantized language, the constraint is 

J drexp{-iq- f}p(f) = | J drexp{-iq- f} £[V x (r)W CT (r)] (27) 

for \q\ < Q and as a shorthand notation we use ^kb,a = also in the following. The proof starts by expressing the 
single-particle operator, 

£(g x V,) exp{-£<f • fj}, (28) 

3," 

in the second-quantized language as 



/ df^ a (r)(qx V)exp{-i(f- r}* CT (r), (29) 



then followed by simple regroupings and the neglect of a surface term. So we find that in the long-wavelength 
approximation we use, the charge density fluctuations exist only if there are vortex excitations in the system. [The 
charge density on the right-hand side (rhs) of Eq. I|27|l is propotional to the vortex densities of the two kinds of fields. 
Strictly speaking the vortex density is defined only by the phase part of a bosonic field, but in the small-momentum 
limit we will neglect the difference. Later a relaxation of the constraint and thereby generation of terms quadratic 
in momenta will be justified by this difference.] If our system is a stable 2D Bose system that would mean that we 
have the case for incompressibility, because to excite a vortex a finite energy is needed so all charge fluctuations are 
supressed. But do we have a stable system? We have to get back to the Hamiltonian expressed in new variables. The 
Hamiltonian is 

W = E |^ - ^ E E PV- exp{-i<f • (fi - mP jt o> + + ^ E A\q)A{q) + V E + V A . (30) 

* =1 ><7 i,<y,],o' \q\<Q \q\<Q 

To eliminate the bare electron mass in the kinetic energy and low energy description we choose, as in Refs. 0, 0], 
that the number of the oscillators is the same as the number of particles, so that the diagonal part (i = j, a = a ') 
of the second term in Eq. I|3l)|l exactly cancels the first kinetic energy term. The bare mass is still present in the 
off-diagonal part of the second term, and if we decompose the q sum as 0, 0] 

£ exp{-iq ■ {n - fj)} = S 2 ^ - fj) - cxp{-iq ■ (f t - fj)}, (31) 
\q\<Q \q\>Q 

we are left with a (5-function interaction among particles and another short-range interaction that may be grouped 
with |J] previously neglected short-range pieces. We are assuming all along that the same kind (layer) of bosons 
(transformed electrons) behave as hard core bosons, so for a — a we see that the S function is ineffective. To 
eliminate the bare mass in the (5-function interaction between the opposite kind bosons we require that they also 
behave mutually as hard core bosons. As we will see, this additional requirement (not due to the fermionic statistics) 
will be very important in the derivation of the low-lying spectrum. 

Therefore, as a result of the transformations made, our Hamiltonian has a free oscillator and Coulomb interaction 
part only. The interaction part in the old variables with the introduced cutoff is 

V = V A + V E = \ P°J d (-q)VA(q)P°J d (q) + l £ p°J d (-q)V E (q)p° M a (q), (32) 

<7,\q\<Q °M<Q 

with Va((1) — 2ire 2 /\q\ and Ve((1) — (27re 2 /|q|) exp{— d\q\}, where d is the distance between layers. If we introduce 

V c {q) = V A (q) + V E (q) and V s (q) = V A (q) - V E (q) , (33) 

and 

P° ld (q) - Pf d (q) + P°i ld (q) and p° ld (q) = pf d (q) - pf(q), (34) 

we can rewrite Eq. I|32() as 

V = E \p° c ld (-<l)VMp c ld (q)+ E \p°s d (-<l)VMp s ld {<i)- (35) 
l<?l<Q ffl<Q 
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To get the expression in new variables for p° ld (q) we should compare Eq. Q24J1. Eq. J2fjJ, and Eq. I|15|l and find 

pf d {q) = q = [A(q)+At(q)] -i^x exp{-i(f ■ r ilCT } (36) 

y/4lT(2s + 1) TT 



and also, as we already found out, we have 



p° s M {q)=p s {q)- (37) 



Further decoupling of the oscillators and particles in V would amount to higher-order corrections to the expressions 
found 0, Q and we can safely neglect the presence of oscillators (terms with A's) in Eq. when discussing the 

low-energy excitations. Then the charge part [the first term in Eq. I|35|l ] can be decomposed into a diagonal and 
off-diagonal part. The diagonal part can be rewritten as 



2m 



i.e., as a kinetic term of particles with mass m c , where 



m, 



E ^ 2 sm 2 % A ., CT (39) 
q\<Q 



Therefore we came to a description of the system in terms of quasiparticles with a mass that is due to interactions. 
These are the expected bosonic dipole objects 0, Q with interaction among them described by the off-diagonal part. 
As overall neutral objects they should make Bose condensate(s) in the ground state and we proceed by taking the 
Bogoliubov expansion of the quasiparticle operators a =|, J, in the second-quantized language as 

*T = a/^ + ^T and *1 = (40) 



where operators r] a ,a =1,1 describe the small fluctuations around the mean field value, where uq is the 

density of particles in each condensate. We also introduce 

Vc = 3L^L and ru = ^ (41) 

new fields that, as we will find out soon, are appropriate for the low-energy description of the system. 
In terms of the new variables, the constraint is, effectively, 

*}* T + - n = n - n + y/nS(nt + Vc) + vile + vhh = 0. (42) 

Please note that the equality here should be understood as the equality of the Fourier transforms of lhs and rhs for q 
small. It is also important to notice that although the constraint effectively is p° ld (q) = for q ^ in the low-energy 
sector and constrains the first term in Eq. (|35J) to vanish, the underlying canonical variables r\ c and r\ s may assume 
nonzero values. To find them, especially r\ s in which we are mostly interested, we do the following decoupling. Due 
to the smallness of ry c and rj s the constraint may be rewritten as 

V f M + Vc(-q)*0 (43) 

so that fields may effectively decouple, satisfying the constraint only approximately. As a result, from the first part in 
Eq. I|35|) . by relaxing the constraint, we get a kinetic term for rj s . There are no other contributions to the second order 
in T] s . From the first part of Eq. (|42|l and in the spirit of the Bogoliubov expansion we may conclude that r\ c is the 
field that couples to the external electromagnetic potential. In our decoupling ansatz r\ c is only very weakly coupled. 
This coincides with the physical picture that we have for bosonic dipoles that (as dipoles) they weakly interact with 
external field and therefore as a system are incompressible Q. 

Applying the Bogoliubov expansion again, and neglecting the difference between hq and n, 

Ps = *}* T - » y/K(nl + Tfc), (44) 
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we are led to the following Hamiltonian for rj s fields, 

n ° = £|^to"H® + j 1 E[v.(-$+vtmv.m.®+Tit(-f)]' (45) 
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As before the hard core boson constraint makes the V s (0) = 2we 2 d part of the interaction ineffective but leaves us with 
7i s that describes an unstable system. Therefore we must impose separately the hard core constraint of composite 
bosons on fields rj s . That amounts locally to the following requirement, 

p2(r) = *j!¥ T + (46) 

where we used the hard boson properties \I/J r \E f cr \E'^\E f cr = "I'J.'J a , f =T an d lj an d ^t^a^-a^-v = 0- Using the 
Bogoliubov expansion, Eq. (|44|) . again this becomes 

n[4(r) +r) s (r)] 2 = n (47) 

that has to be imposed on rj s fields. Note that here we also used the incompressibility property of the system in 
the low-energy region for the charge degrees of freedom, on the rhs irrespective of the length scale [Eq. (|42|l Fourier 
transformed for any q\ . We had to make this assumption because we are incorporating a piece of short-range physics 
into the long-wavelength description. Please also note that this is an operator identity, where the automatic neglect 
of the quadratic terms on the lhs of the equation, in the Bogoliubov expansion, is not allowed. 

The constraint we handle in the usual way, switching to the Lagrangian formulation with fields r) s , 77] and a field A 
that enforces the constraint [TT| . The generating functional is 

Z = J V Vs J VritJ VX expj-^rfr J d 2 x (4d T r, s + H s {x) + {[ Va (x) + 4{x)] 2 - 1} iX(x, r))| , (48) 

where 

H s (x,t) = ^-fvt ■ V77, + 7 / d 2 y[i ls {x) + vt(x)]V s (x - ff)[r] s (y) + 4{y)]. (49) 
2m c 4 J 

The constraint approximately commutes with the Hamiltonian in the long-wavelength limit (using this property we 
combined contributions into a single exponential), and so we will take A, r (imaginary time) to be independent. Also, 
at the mean- field level, we are allowed to assume that A is space independent. 
Introducing Bogoliubov transformations on rj s (q,T) fields, 

rj s (q, t) = ot(q) exp{iw 9 T} cosh^+ cr(— q) exp{— iuj q r} sinh^ (50) 

[where ct(q) and c^{q) are new canonical fields], we get after standard transformations that diagonalize the problem 
(see also Ref. [l2|); the following mean-field expression for Z: 



where 



rfA n T3^^ exp { -^ (A)}, (51) 



e(q, A) = ^) 2 + K( 9 ) + 4 8 A]e§, (52) 
with e c q — q 2 /2m c , and the domain of <f s is again the disk with radius Q. Also 

£o(A) = i^[e( g ,A)-4]-iA£ (53) 
9 g 

with the q summations where the cutoff Q is understood. The e(q, A)'s are the usual Bogoliubov energies, the results 
of the Bogoliubov transformation, now requiring also a suitable A to get the final expression for the mode dispersion 
we are looking for. We approximate A in the saddle point approach (see, for example, Ref. searching for a 

stationary point of F(T, A), from the following expression for Z m f , 



Zmf = J rfAexp{-/3F(T, A)}, 



(54) 
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FIG. 1: The stationary point values for the velocity of the out-of-phase mode (o), the values from the pseudospin theory in the 
Hartree-Fock approximation from Ref. (A), and exact diagonalization results (*) from Ref. ||. 



i.e., look for the solution of 



which effectively becomes 



dF(T, A) 
dX 



dE Q (X) 
dX 



= 0, 



0, 



(55) 



(56) 



in the T — > limit. 

We solved the equation numerically finding only solutions Ao with i Xq real and positive (therefore, as usual |13|. we 
found a path and a saddle point in the complex A plane), and results are depicted in Fig. 1. In Fig. 1 we also plotted 
the Hartree-Fock result of Ref. and the same reference exact diagonalization results at d = 0.51b and d = lis- 
To a good approximation we can claim a linear dependence for small d of the Bogoliubov velocity though with the 
values significantly reduced from the Hartree-Fock results. But at d = 0.51b all three data points are very close to 
each other. For larger values of d, d ~ 11b and larger, both approximation schemes fail to capture the quantum 
fluctuations that increase with d |8(- 

The result for e(q), the dispersion of the out-of-phase mode in the d = case, 



e(<7) = e * = w 

Hi and 0- Namely e(q) 



(57) 



e l ls\V 7 iq 2 while the Hartree-Fock result is 
eHni^pK /8)y/2q 2 . It is interesting to note that if we use the expression for m c conjectured in the generalized theory 



is fairly close to the estimate of Refs. 

0,13 that includes higher-momentum physics, we exactly get the Hartree-Fock result. 

It is also interesting to speculate about the discrepancy between our and the Hartree-Fock result for small d. Part of 
it might be due to our low- energy, low-momentum limited approach, but it might also well be due to the incompletness 
of the underlying analogy 0, 1 1 -ll ] of the d ^ system description compared to the one of a repulsively interacting 
Bose gas. (The analogy of the d = case to an ideal Bose gas is complete, as we found out.) The incompletness 
might follow from the modifications of the composite boson picture due to the presence of composite fermions as 
proposed in Ref. |l5j |. The composite fermions come into relevance very soon as d acquires a nonzero value, and their 
number rapidly increases with d |l5j | . If we are allowed to view bosons and fermions to a first approximation as weakly 
interacting through a short-range interaction (more precisely here interacting are differences between up and down 
bosons and fermions, respectively) we can borrow considerations applied to the Bose-Fermi mixtures in optical traps; 
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see, for example, Refs. and [l7|. When we have a small number of fermions the Bogoliubov mode can, in fact, 
enhance its value (on the other hand, decrease) and damping follows when fermions proliferate. In that sense we can 
expect that by considering fluctuations around our mean field solution (which maps the problem to an interacting 
bose gas with a Bogoliubov mode), the value of the Bogoliubov velocity may increase for small d, stay almost the 
same for intermediate values, and decrease and even acquire damping for large d. In this sense here we have set up a 
formalism necessary to check these considerations in a future work. 

It is important to address the case when we do not have the same density of particles in the up and down layer 
0,0]. Then in general we have, instead of Eq. 

P°s d =Ps~ Ot - x ex P{-*9 - (58) 

when, as before, we neglected the magnetoplasmon part. [ and v\ (y^ + = 1) are the filling factors of each layer 
separately] For p a , a =f, j in p s = p^ — p^ we assume the following form: 

Ax = E^ x exp{-i<f- r^ a ] + C ^(q- Pj,a) exp{-i<f • (59) 

j j 

The second part is the longitudinal component of the paramagnetic current, and the term should appear, in general, 
when compressible low- lying degrees of freedom are present. If C in Eq. I|59|l is the same for both layers then 
p c = P]- + Pi, because the total component of the current is zero in the charge channel due to its incompressible 
nature. Substituting Eq. (|58|l with Eq. (|59|) in the projected Hamiltonian [Eq. (|35|l ] and collecting all diagonal terms 
of the dipole expansion, we get 

J- = £ ,Vfe,J - ^(^A \ (- + -), (60) 

m Su ^ Q I 2 2 2 J 4 V^T v i) 

as a generalization of Eq. I|39|) to the case of imbalanced layers. This is the mass of the r\ s field defined as r\ s — 
\Jv\ r \X ~ \ / ^l r ll- So we assumed that we can apply the Bogoliubov theory and with neglect of some residual dipole- 
dipole interaction in the pseudospin channel, our problem reduces to the one expressed in Eq. (|i*S|> and Eq. l|4^|) 
in which instead of the mass, ro c , we have msv The assumption is based on the expectation that the pseudospin 
channel is compressible. For not large — fjj = 8v, the velocity of the Bogoliubov mode decreases quadratically 
with 5v as a consequence of Eq. (ffj()|) . in agreement with Ref. 20]. A more detailed investigation of the influence of 
the dipole-dipole interaction is needed for general 6v. 

We would like to address also the case of huge imbalance, when we take, for example, 3> vi |2l|. As — > 
0, l/mg u — > oo, and r/ s = ^Jv^r\\ — y/v\f)i ~ Vh which, probably signals the incompressible physics of the | layer. 
Therefore, to find out more about the physics of J, quasiparticles, we must go back to the beginning formulation, and 
apply a different decomposition. Namely we will take ( in the limit i/| > j/j) 

p oid w p oid _ 2pi = -^(qx p^)exp{-iq- f^} - 2 Pl . (61) 

It is easy to see that the first term in Eq. I|61|) would lead to the effective mass for all quasiparticles, in the first 
approximation, equivalent to what we would have if there was only one single layer with v = 1. Next considering the 
cross term, 

^2[-2 Pi (- g )]^p° w ( (Z ), (62) 

q<Q 

and taking the expression in Eq. (|59|l for we get from Eq. I|62|) for the J, quasiparticle mass, 

«E^Vsin 2 % p ^. (63) 

me <3<Q 

Because of the incompressible | background we can neglect | and j cross terms. If we again also assume irrele- 
vance of the remaining | dipole-dipole interactions (for the low- momentum physics), our effective Hamiltonian for j 
quasiparticles is 

^Ej^t^ + E*^*- (64) 
q<Q e q<Q 
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In Ref. |2lJ a physical picture of a Bose gas of excitons and dipoles with density n± was developed for the case 
v *\ ^ v i- The mass m e we derived is the result of the low-momentum theory. As in the d = 0, = case we expect 
that in the generalized theory 0, El the cutoff in Eq. would be replaced by a Gaussian in the momentum space 
and the expression would coincide with the one in Ref. |2l|. In this sense with the assumption made, also in this case 
due to the comparison to Ref. l2l[. we can claim a complete analogy to a weakly interacting Bose gas in d — ► limit. 

References [y, 0, |22] and Ref. [231 (when not considering spiral states) agree on the dependence of the Bogoliubov 
velocity (in the n-f = case). Possible additions of quantum fluctuations to this value can be extracted from Ref. 
|24| . There, due to the justified assumption of the suppression of charge fluctuations; a Schwinger boson mean field 
theory was used with the requirement on the single occupancy of the Schwinger boson in a lowest-Landau-level basis. 
In this work we were primarily concerned with the establishment of the concept of a composite boson, and we only 
set up the stage for considering fluctuations beyond generalized CS mean-field theory that is based on this concept. 
(A composite boson approach may prove useful for the study of quantum phase transitions in the bilayer |l5l l25| , and 
building of the physical picture of the bilayer in analogy with the picture based on composite fermions in the single 
layer.) Our mean-field theory and the usual theory do not agree somewhat, although they agree to a much better 
degree than the usual CS theory 0,113 (linear dependence on small d and the absence of the bare mass). Inclusion of 
the fluctuations in our hard-core (belonging to different layer) CS boson model that is probably related to the model 
with the single occupancy of Schwinger bosons and comparison to Ref. |24j ] are planned for future work. It would 
be important to probe the significance of the fluctuations around d ~ Ib- Any strong instability of the Bogoliubov 
mode velocity would signal, in the composite-boson- composite- fermion model (see above and Ref. |15|L the phase 
separation of the two fluids [23] and the proposed first-order transition [2^,|23. Then, from the composite-boson 
point of view, we would be able to address in more detail the extraodinary experiments done on the bilayer |30l l3ll ] . 
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